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Abstract 

In this paper, we consider a cognitive radio (CR) network in which the unlicensed (secondary) 
users are allowed to concurrently access the spectrum allocated to the licensed (primary) users provided 
that their interference to the primary users (PUs) satisfies certain constraints. We study a weighted 
sum rate maximization problem for the secondary user (SU) multiple input multiple output (MIMO) 
broadcast channel (BC), in which the SUs have not only the sum power constraint but also interference 
constraints. We first transform this multi-constraint maximization problem into its equivalent form, which 
involves a single constraint with multiple auxiliary variables. Fixing these multiple auxiliary variables, 
we propose a duality result for the equivalent problem. Our duality result can solve the optimization 
problem for MIMO-BC with multiple linear constraints, and thus can be viewed as an extension of the 
conventional results, which rely crucially on a single sum power constraint. Furthermore, we develop 
an efficient sub-gradient based iterative algorithm to solve the equivalent problem and show that the 
developed algorithm converges to a globally optimal solution. Simulation results are further provided 
to corroborate the effectiveness of the proposed algorithm. 
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I. Introduction 

Cognitive radio (CR), as a promising technology to advocate efficient use of radio spectrum, 
has been a topic of increasing research interest in recent years [l]-[7]. CR allows an unlicensed 
(secondary) user to opportunistically or concurrently access the spectrum initially allocated to the 
licensed (primary) users provided that certain prescribed constraints are satisfied, thus having 
a potential to improve spectral utilization efficiency. In this paper, we study a weighted sum 
rate maximization problem for the secondary user (SU) multiple input multiple output (MIMO) 
broadcast channel (BC) in a concurrent CR network, in which the SUs have not only the sum 
power constraint but also interference constraints. 

A. System Model and Problem Formulation 

With reference to Fig. 1, we consider the A'-SU MIMO-BC with Nt transmit antennas and 
Nr receive antennas in a CR network, where the K SUs share the same spectrum with a single 
primary user (PU) equipped with one transmitter and one receiver^ The transmit-receive signal 
model from the BS to the zth SU denoted by SUj, for i = 1, . . . ,K, can be expressed as 

y, = HiX + Zi, (1) 

where is the A^^ x 1 received signal vector. Hi is the Nr x Nt channel matrix from the BS to 
the SUj, X is the Nt x 1 transmitted signal vector, and Zi is the Nr x 1 Gaussian noise vector 
with entries being independent identically distributed random variables (RVs) with mean zero 
and variance cx^. Consider ho as the A'^t x 1 channel gain vector between the transmitters of the 
BS and the PU. We further assume that for i = 1, . . . , K, and ho remain constant during 
a transmission block and change independently from block to block, and Hi for i = 1, . . . , K, 
and ho are perfectly known to the BS and SUj. This requires that the SUs can "cognitively" 
obtain the information of its neighboring environment. In practice, certain cooperation in terms 
of parameter feedback between the PU and the BS is needed. To achieve that, the protocol for 
the SU network can be designed as follows: every frame contains sensing sub-frame and data 
transmission sub-frame. During the sensing sub-frame, BS can transmit training sequences to 
SUs as well as to the PU so that the SUs can estimate the channel matrix Hi, and the PU can 
measure the vector ho- After that, this information will be sent back to the BS via a feedback 
channel. 

'Expect for explicitly stated, we restrict our attention to a single PU case in the rest of this paper for convenience of description. 
The results derived for the single PU case can be readily extended to the multiple PU case, which is discussed in Remark 4. 
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We next consider the weighted sum rate maximization problem for the Ji-SU MIMO-BC 
in a CR network, simply called the CR MIMO-BC sum rate maximization problem, which, 
mathematically, can be formulated as 

Problem 1 (Main Problem): 

K 

max Wir° (2) 

K K 

subject to ^hlQ\ho < Pt, and ^tr(Q^) < P„, 

1=1 i=l 

where is the rate achieved by SUj, Wi is the weight of SU/, denotes the Nt x Nt transmit 
signal covariance matrix for SUj, h denotes that is a semidefinite matrix, Pt denotes the 
interference threshold of the PU, and P„ denotes the sum power constraint at the BS. In a non- 
CR setting, similar weighted sum rate optimization problems for the multiple input single output 
(MISO) BC and the MIMO-BC have been studied in [8] [9], respectively. The key difference is 
that in addition to the sum power constraint, an interference constraint is applied to the SUs in 
the CR MIMO-BC, i.e., the total received interference power 

E£i hlQ^iho at the PU is below 

the threshold Pt. 

Remark 1: It has long been observed that the optimal sum rate for MIMO BC with a single 
sum power constraint is equal to the optimal sum rate of the dual MIMO multiple access channel 
(MAC) with the same sum power constraint [10] [11] [12]. However, this conventional BC-MAC 
duality can only be applied to the case with a single sum power constraint (even not applicable 
to an arbitrary linear power constraint). Hence, the additional interference power constraint in 
Problem 1 makes the existing duality cannot be applied. The new duality result proposed in this 
paper generalizes the previous results as special cases. Moreover, it is worth to note that any 
boundary point of the capacity regions of the MIMO-MAC and the MIMO-BC can be expressed 
as a weighted sum rate for a certain choice of weights [13] [14]. Thus, by varying the weights 
of the SUs in Problem 1, the entire capacity region of the CR MIMO-BC can be obtained. 

B. Related Work 

The present paper is motivated by the previous work on the information-theoretic study of the 
MIMO-BC under a non-CR setting. It has been shown in [11] [12] [15] that under a single sum 
power constraint, the sum-capacity of the non-CR MIMO-BC can be achieved by the dirty paper 
coding (DPC) scheme. Furthermore, the paper [16] shows that the rate region achieved by the 
DPC scheme is indeed the capacity region of such a channel. However, the power allocation and 
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beamforming strategies to achieve the capacity region have been not considered in these papers. 
Moreover, it has been shown in [17] [18] that under the single sum power constraint, the equally 
weighted sum rate maximization problem, simply called the sum rate problem, for the MIMO-BC 
can be solved by solving its dual MIMO MAC sum rate problem, which is also subject to a single 
sum power constraint. In [17], a cyclic coordinate ascent algorithm was proposed to solve the dual 
MIMO-MAC problem while in [18] this sum-power constrained dual problem was decoupled 
into an individual-power constrained problem, which can be solved by using an iterative water- 
filling algorithm [19]. Even though these algorithms proposed in [17] [18] can solve the sum 
rate optimization problem for the non-CR MIMO-BC via the MAC-BC duality, they are not 
applicable to the general weighted sum rate problem. In [8], a generalized iterative water-filling 
was proposed to solve the weighted sum rate problem for the MISO-BC where each user has a 
single receive antenna. However, the proposed algorithm is not applicable to the general MIMO- 
BC case. Furthermore, an efficient algorithm was proposed to solve the MIMO-BC weighted 
sum rate problem with a single sum power constraint in [9]. These aforementioned results are 
based on the conventional BC-MAC duality, which cannot be applied to solve the weighted sum 
rate problem with multiple constraints (the case of interest in this paper). Recently, the paper [20] 
investigated a different MIMO-BC weighted sum rate maximization problem which is subject 
to per-antenna power constraints instead of the single sum power constraint, and established a 
new minimax duality which is different from the conventional BC-MAC duality. A Newton's 
method based algorithm was proposed to solve this minimax problem. In this paper, we consider 
a more general case where the power is subject to multiple linear constraints instead of the sum 
power constraint or per-antenna power constraints, and propose a new BC-MAC duality result 
to extend the conventional duality result so that it can solve the problem with multiple arbitrary 
linear constraints. A Karush-Kuhn-Tucker (KKT) condition based algorithm is developed to solve 
the problem. 

C. Contribution 

Throughout the paper, we consider the CR MIMO-BC weighted sum rate maximization 
problem as defined in Problem 1. As the main contribution of this paper, our solution is 
summarized in the following. 
1) We prove that in the CR MIMO-BC, the multi-constraint weighted sum rate maximization 
problem (Problem 1) is equivalent to a single-constraint weighted sum rate maximization 
problem with multiple auxiliary variables. 
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2) For the equivalent problem, we establish a duality between the MIMO-BC and a dual 
MIMO-MAC when the multiple auxiliary variables are fixed as constant. This duality is 
applicable to MIMO-BC with arbitrary linear power constraint, and can be viewed as an 
extension of the conventional MIMO MAC-BC duality result [10] [11] [12], which is only 
valid for the problem with a single sum power constraint. 

3) For the weighted sum rate maximization problem of the dual MEMO MAC, the existing 
iterative water-filling based algorithm [17], [18] is not applicable. We propose a new primal 
dual method based iterative algorithm [21] to solve it. Furthermore, we propose a sub- 
gradient based iterative algorithm to solve the main problem of the paper. Problem 1, and 
show that the proposed algorithm converges to the globally optimal solution. 

D. Organization and Notation 

The rest of the paper is organized as follows. In Section II, we transform the CR MIMO-BC 
weighted sum rate maximization problem (Problem 1) into its equivalent form, and introduce a 
MAC-BC duality between a MIMO-BC and a dual MIMO-MAC. Section in presents an primal 
dual method based iterative algorithm to solve the dual MIMO-MAC weighted sum rate problem. 
In Section IV, a MAC-BC covariance matrix mapping algorithm is proposed. Section V presents 
the complete algorithm to solve the CR MIMO-BC weighted sum rate maximization problem. 
Section VI provides several simulation examples. Finally, Section VII concludes the paper. 

The following notations are used in this paper. The boldface is used to denote matrices and 
vectors, and (■)^ denote the conjugate transpose and transpose, respectively; Im denotes an 
M X M identity matrix; tr(-) denotes the trace of a matrix, and \x\^ denotes max(x, 0); (■)'' and 
(■)™ denote the quantities associated with a broadcast channel and a multiple access channel, 
respectively; E\-] denotes the expectation operator. 

II. Equivalence and Duality 

Evidently, the MIMO-BC weighted sum rate maximization problem under either a non-CR 
or a CR setting is a non-convex optimization problem and is difficult to solve directly. Under a 
single sum power constraint, the weighted sum rate problem for MIMO BC can be transformed 
to its dual MIMO MAC problem, which is convex and can be solved in an efficient manner 
[8] [9]. In the CR setting, the problem (Problem 1) has not only a sum power constraint but 
also an interference constraint. The imposed multiple constraints render difficulty to formulate 
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an efficiently solvable dual problem. To overcome the difficulty, we first transform this multi- 
constrained weighted sum rate problem (Problem 1) into its equivalent problem which has a 
single constraint with multiple auxiliary variables, and next develop a duality between a MIMO- 
BC and a dual MIMO-MAC in the case where the multiple auxiliary variables are fixed. 

A. An Equivalent MIMO-BC Weighted Sum Rate Problem 

In the following proposition, we present an equivalent form of Problem 1 (see Appendix B 
for the proof). 

Proposition 1: Problem 1 shares the same optimal solution with 
Problem 2 (Equivalent Problem): 

K 

min max Wir\ (3) 

K K 

subject to qt{Y, ^WiK - Pt) + QuiY, ^^"(Q') - ^ 0' ^ 
1=1 1=1 
where qt and g„ are the auxiliary dual variables for the respective interference constraint and 

sum power constraint. 

It can be readily concluded from the proposition that the optimal solution to Problem 2 also 
satisfies XliLi ^iQi < Pt and J2iLi ) < Pu simultaneously since it is also the optimal 
solution to Problem 1 . Finding an efficiently solvable dual problem for Problem 2 directly is still 
difficult. However, as we show later, when qt and g„ are fixed as constants. Problem 2 reduces 
to a simplified form, which we can solve by applying the following duality result. 

B. CR MIMO BC-MAC Duality 

For fixed qt and g^. Problem 2 reduces to the following form 
Problem 3 (CR MIMO-BC): 

K 

max WiT-^i (5) 

K K 

subject to qt hlQ\K + quY, < p. (6) 

i=l i=l 

where P := qtPt + quPu- Since qt and g„ are fixed, P is a constant in Problem 3. The constraint 
(6) is not a single sum power constraint, and thus the duality result established in [17] is not 
applicable to Problem 3. Therefore, we formulate the following new dual MAC problem. 
Proposition 2: The dual MAC problem of Problem 3 is 
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Problem 4 (CR MIMO-MAC): 

K 

max > Wirf (7) 

{gr}f=i: Q™bo ^ 

subject to ^tr(gf)a2<P, (8) 

i=l 

where is the rate achieved by the ith user of the dual MAC, Q™ is the transmit signal 
covariance matrix of the zth user, and the noise covariance at the BS is qthoh^ + QulNf 

Remark 2: According to Proposition 2, for fixed qt and Qu, the optimal weighted sum rate of 
the dual MAC is equal to the optimal weighted sum rate of the primal BC. From the formulation 
perspective, this duality result is quite similar to the conventional duality in [10] [11] [12]. 
However, as shown in Fig. 2, one thing needs to highlight is that the noise covariance matrix 
of the dual MAC is a function of the auxiliary variable qt and g^, instead of the identity matrix 
[12]. This difference comes from the constraint (6), which is not a sum power constraint as in 
[12]. Note that when qt = 0, the duality result reduces to the conventional BC-MAC duality in 
[12]. 

As illustrated in Fig. 2, Proposition 2 describes a weighted sum rate maximization problem for 
a dual MIMO-MAC. To prove the proposition, we first examine the relation between the signal 
to interference plus noise ratio (SINR) regions of the MIMO-BC and the dual MIMO-MAC. 
Based on this relation, we will show that the achievable rate regions of the MIMO-BC and the 
dual MIMO-MAC are the same. 

In the sequel, we first describe the definition of the SINR for the MIMO-BC. It has been 
shown in [16] that the DPC is a capacity achieving scheme. Each set of the transmit covariance 
matrix determined by DPC scheme defines a set of transmit and receive beamforming vectors, 
and each pair of these transmit and receive beamforming vectors forms a data stream. In a 
beamforming perspective, the BS transmitter have Nt x K beamformers, ■Uj j, for i = 1, ■ ■ ■ , ii', 
and j = 1, ■ ■ ■ , A^j. Therefore, the transmit signal can be represented as 

K Nt 

X = ^ ] ^ ] ^iyj'^iyjl 
i = l j = l 

where Xi^j IS 3. scalar representing the data stream transmitted in this beamformer, and E[xf ] = 
Pi J denotes the power allocated to this beamformer. At SUj, the receive beamformer correspond- 
ing to Uij is denoted by Vij. The transmit beamformer Uij and the power pij can be obtained 
via the eigenvalue decomposition of Q^, i.e., = U\PiUi, where Ui is a unitary matrix, and 
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Pi is a diagonal matrix. The transmit beamformer Uij is the jth column of C/j, and is the 
jth diagonal entry of Pi. With these notations, we express the SINR^^ as 

, Pi Au] -H^Vi A"^ 
SINRl", = ^ , ' , , . (9) 

It can be observed from (9) that the DPC scheme is applied. This can be interpreted as follows. 
The signal from SUi is first encoded with the signals from other SUs being treated as interference. 
The signal from SU2 is next encoded by using the DPC scheme. Signals from the other SUs will 
be encoded sequentially in a similar manner. For the data streams within SU^, the data stream 1 
is also encoded first while the other data streams are treated as the interference. The data stream 
2 is encoded next. In a similar manner, the other data streams will be sequentially encoded. The 
encoding order is assumed to be arbitrary at this moment, and the optimal encoding order of 
Problem 2 will be discussed in Section III. 

To explore the relation of the SINR regions of the dual MAC and the BC, we formulate a 
following optimization problem 

K K 

min qt V hlQ\K + g« V iv{Q\) - P 

{Q!}f.,:Q^bo tr tr (10) 

subject to SINR- ^ > 7^ ,,, 
where 7^ denotes the SINR threshold of the jth data stream within the SU^ for the BC. Note 
that the objective function in (10) is a function of signal covariance matrices and the constraints 
are SINR constraints for the K-S\J MIMO-BC. 

It has been shown in [20] and [22] that the non-convex BC sum power minimization problem 
under the SINR constraints can be solved efficiently via its dual MAC problem, which is a 
convex optimization problem. By following a similar line of thinking, the problem in (10) can 
be efficiently solved via its dual MAC problem. Similar to the primal MIMO-BC, the dual 
MIMO-MAC depicted in Fig. 2 consists of K users each with A'^,. transmit antennas, and one 
BS with Nt receive antennas. By transposing the channel matrix and interchanging the input and 
output signals, we obtain the dual MIMO-MAC from the primal MIMO-BC. For the covariance 
matrices of the dual MIMO-MAC, we apply the eigenvalue decomposition, 

Nr 

Ql = ViX^^\ = Y,%3^^A,p (11) 

where is the jth column of Vi, and qij is the jth diagonal entry of Aj. For user i, Vij 
is the transmit beamforming vector of the jth data stream, the power allocated to the jth data 
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Stream equals qij, and the receive beamforming vector of the jth data stream at the BS is Uij. 
The SINK of the dual MIMO-MAC is given by 



SINK- = . . , ^ q^MJH^v.,J^ ^ ^^2) 

where Rw := qtRo + Qui Nt is the noise covariance matrix of the MIMO-MAC with Ro := hoh\. 
In the dual MIMO-MAC, R^ depends on qt and defined in (10) whereas the noise covariance 
matrix in the primal MIMO-BC is an identity matrix. It can be observed from (12) that the 
successive interference cancelation (SIC) scheme is used in this dual MIMO-MAC, and the 
decoding order is the reverse encoding order of the primal BC. The signal from SUi^ is first 
decoded with the signals from other users being treated as interference. After decoded at the 
BS, the signals from SU^: will be subtracted from the received signal. The signal from SUi^_i 
is next decoded, and so on. Again, the data streams within a SU can be decoded in a sequential 
manner. 

For the dual MIMO-MAC, we consider the following minimization problem similar to the 
problem (10) 

K 

min Vtrfgn^^-^ 
{QTlfir Qrbo ^ (13) 

subject to SINR^j > 7,^. 
The following proposition describes the relation between the problems (10) and (13). 

Proposition 3: For fixed qt and g„, the MIMO-MAC problem (13) is dual to the MIMO-BC 
problem (10). 

Proof: The constraints in (10) can be rewritten as 

I t trt 12 K Nr Nr 

"•■''"■f^-""' > E Y.P''M,Hlv.f+ E PM,Hlv.,f + a^. (14) 

^'^'^ k=i+l 1=1 l=j+l 

Thus, the Lagrangian function of the problem (10) is 



Li{Q\, . . . , Q^, A 



K K K Nr I t trt 12 

-q, Y: hlQ\K + E tr(Q^) - P - E E 

i=l i=l i=l j=l "'-^ 

K Nr Nr 

- E Y.p''M,iHiv,,\' - E p^,HlH^v.J' - (15) 

k=i+l 1=1 l=j+l 
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K Nr K Nr \ trt t u 



i-l Nr j-l 



1: 



- X] K,iH\vk,iv\ iHk - ^ \iH\vijv\iHi - R^j Uij, (16) 

k=l 1=1 1=1 

where Aj j is the Lagrangian mukiplier. Eq. (16) is obtained by applying the eigenvalue decom- 
position to and rearranging the terms in (15). The optimal objective value of (10) is 

max min Li{Q\, . . . ,Q)^, Xij). (17) 
Q\,-,Qk 

On the other hand, the Lagrangian function of the problem (13) is 

K K Nr Jji t JJ 

L^iQT, . . . , Q^, s.„) = j:y. - ^ - E E kM /'^' 'I;:' ' 

1=1 j=l i=l j=l ^^'^ 

i-l Nr j-l 

1k,iHlvk,ivl iHk - E %iH\vi^iv\iHi - R^)uij, (18) 

k=l 1=1 1=1 

where 5i j is the Lagrangian multiplier. Eq. (18) is also obtained by applying eigenvalue decom- 
position to Q™. The optimal objective value of (13) is 

max min L2(Q^ . . . , Q^, 5„). (19) 

Note that if we choose qij = Xij, Sij = pij, and the same beamforming vectors Uij and Vij 
for both problems, (16) and (18) become identical. This means that the optimal solutions of (17) 
and (19) are the same. ■ 

Proposition 3 implies that under the SINR constraints, the problems (10) and (13) can achieve 
the same objective value, which is a function of the transmit signal covariance matrices. On the 
other hand, under the corresponding constraints on the signal covariance matrix, the achievable 
SINR regions of the MIMO-BC and its dual MIMO-MAC are the same. Mathematically, we 
define the respective achievable SINR regions for the primal MIMO-BC and the dual MIMO- 
MAC as follows. 

Definition 1: A SINR vector 7 = (71,1, . . . , 7i,Art, . . . , 7i^,Art) is said to be achievable for 
the primal BC if and only if there exists a set of Q'l, . . . , Q\ such that qt Y^=i ^^IQ^^^o + 
Qu J2f=i ) — P < C for a constant C and the corresponding SINR^^ > 7^^. An achievable 
BC SINR region denoted by TZbc^ is a set containing all the BC achievable 7. 

Definition 2: A SINR vector 7 = (71,1, . . . , 7i,Art, • • • , lK,Nt) is said to be achievable for the 
dual MAC if and only if there exists a set of , . . . , such that Ylf=i ^AQT)<^^ - P < C for 
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a constant C and the corresponding SINR™ > 7^ An achievable MAC SINR region denoted 
by TZmac, is a set containing all the MAC achievable 7. 
In the following corollary, we will show TZmac = T^bc- 

Corollary 1: For fixed qt and q.^, and a constant C, the MIMO-BC under the constraint 
E£i hlQ^iho + g„ E£i tr(Qi ) - P < C and the dual MIMO-MAC under the constraint 
E£i tr(Qr)t^^ -P <C achieve the same SINR region. 

Proof: For any 7 G 7^mac> by Definition 2, there exists a set of Q^, . . . , such that 
Ei=i ^AQTW - P < C and the corresponding SINR^^. > 7^^. It can be readily concluded from 
Proposition 3 that there exists a set of QJ, ... , Q^^- such that qt Ylf=i h\Q\ho^qu Yl!i=i ^^(Qi ) — 
P < C and the corresponding SINR^j > 7,j This implies 7 G T^bc- Since 7 is an arbitrary 
element in T^mac, we have TZmac ^ T^bc- In a similar manner, we have TZ^c ^ T^mac- The 
proof follows. ■ 
We are now in the position to prove Proposition 2. 

Proof of Proposition 2: According to Corollary 1, if C = 0, then under the constraint 
Qt Eti hlQlK + qu Ef=i tr(Q' ) < P for the BC and the constraint J2f=i <Q7W < P 
for the dual MAC, the two channels have the same SINR region. Since the achievable rates 
of user i in the MIMO-MAC and the MIMO-BC are = Ef=ilog(l + SINR™) and r^^ = 
X]j^i log(l + SINR^j), the rate regions of the two channels are the same. Therefore, Proposition 
2 follows. ■ 

Note that due to the additional interference constraint. Problem 2 cannot be solved by using 
the established duality result in [11] and [12], in which only a single sum power constraint was 
considered. Our duality result in Proposition 2 can be thought as an extension of the duality 
results in [11] [12] to a multiple linear constraint case. Moreover, as will be shown in the 
following section, our duality result formulates a MIMO-MAC problem (Problem 4), which can 
be efficiently solved. 

III. Dual MAC Weighted Sum Rate Maximization Problem 

In this section, we propose an efficient algorithm to solve Problem 4. With the SIC scheme, 
the achievable rate of the kth user in the dual MIMO-MAC is given by 

Tk = log fell In t • ^^^^ 

For the MIMO-MAC, the equally weighted sum rate maximization is irrespective of the decoding 
order. However, in general the weighted sum rate maximization in the MIMO-MAC is affected 
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by the decoding order. We thus need to consider the optimal decoding order of the SIC for the 
dual MIMO-MAC, and further need to consider the corresponding optimal encoding order of 
the DPC for the primal BC. 

Let TT be the optimal decoding order, which is a permutation on the SU index set {1, ■ ■ ■ , K}. 
It follows from [14] that the optimal user decoding order vr for Problem 4 is the order such that 
Wtt{i) > ^7r(2) > ■ ■ ■ > is satisfied. The following lemma presents the optimal decoding 

order of the SIC for the data streams within a SU (see Appendix C for the proof). 

Lemma 1: The optimal data stream decoding order for a particular SU is arbitrary. 

Due to the duality between the MIMO-BC and the MIMO-MAC, for Problem 3, the optimal 
encoding order for the DPC is the reverse of vr. Because of the arbitrary encoding order for the 
data streams within a SU, if we choose a different encoding order for the BC, the MAC-to-BC 
mapping algorithm can give different results which yield the same objective value. Hence, the 
matrix achieving the optimal objective value are not unique. With no loss of generality, we 
assume wi > W2 > ■ ■ ■ > wk for notational convenience. 

According to (20), the objective function of Problem 4 can be rewritten as 

K i 

fiQi, ■ ■ ■ , Q^) := 5^ A, log \R^ + J2 HjQfH% (21) 

i=i j=i 

where Aj := Wi — Wi+i, and wk+i := 0. Clearly, Problem 4 is a convex problem, which can be 
solved through standard convex optimization software packages directly. However, the standard 
convex optimization software does not exploit the special structure of the problem, and thus is 
computationally expensive. An efficient algorithm was developed to solve a weighted sum rate 
maximization problem for the SIMO-MAC in [8]. However, since this algorithm just consider 
the case where each users has a single data stream, it is not applicable to our problem. In the 
following, we develop a primal dual method based algorithm [21] to solve this problem. 
We next rewrite Problem 4 as 

K 

mRx f{Q^, ) subject to J] tr{QT) < P. (22) 

{Qi }i=i- Qi ^—i 

Recall that the positive semi-definiteness of is equivalent to the positiveness of the eigen- 
values of Qf, i.e., Qij > 0. Correspondingly, the Lagrangian function is 

K K Mi 

L{Q^, ■■■,Ql,X, 5.,) = f{Q^, . . . , ) - A ( ^ tr{QT) - P) ^^'^^^^^^ (^3) 

i=l i=l j=l 
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where A and 5ij are Lagrangian multipliers. According to the KKT conditions of (22), we have 



dQT 



-XIN. + Y,klE^,J=0, (24) 

K 

\{Y,HQT) - P) = 0, (25) 



i=l 

kjlij = 0' (26) 

where Eij := dqij/dQ^. Notice that it is not necessary to compute the actual value of 5ij and 
Eij, because if 5ij ^ 0, then = 0. Thus, the semi-definite constraint turns into qij = [qij]'^- 
Thus, we can assume 5ij = 0. 

The dual objective function of (22) is 

(7(A) = ^ max L(Q- ■■■,Q-,A). (27) 

Because the problem (22) is convex, it is equivalent to the following minimization problem 

min g{X) subject to A > 0. (28) 

A 

We outline the algorithm to solve the problem (28). We choose an initial A and compute the 
value of g{\) (27), and then update A according to the descent direction of g{\). The process 
repeats until the algorithm converges. 

It is easy to observe that all the users share the same A, and thus A can be viewed as a water 
level in the water filling principle. Once A is fixed, the unique optimal set {Q™, . . . , Q'^} can be 
obtained through the gradient ascent algorithm. In each iterative step, Q'^ is updated sequentially 
according to its gradient direction of (23). Denote by Qf{n) the matrix at the nth iteration 
step. The gradient of each step is determined by 

"^QT^-- 9Qr(n - 1) -^^^^ ^^^^ 

Thus, Qf{n) can be updated according to 

QTin) = [QTin - 1) + tvglL 

where t is the step size, and the notation is defined as [A]+ := with Aj and 

Vj being the jth eigenvalue and the corresponding eigenvector of A respectively. The gradient 
in (29) can be readily computed as 



3=k 



(30) 
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where Fj{Q'^, ■ ■ ■ , Q^) := R.^^ + Y^l^i HlQfHi. We next need to determine the optimal A. 
Since the Lagrangian function g(\) is convex over A, the optimal A can be obtained through 
the one-dimensional search. However, because g(\) is not necessarily differentiable, the gradient 
algorithm cannot be applied. Alternatively, the subgradient method can be used to find the optimal 
solution. In each iterative step, A is updated according to the subgradient direction. 

Lemma 2: The sub-gradient of g{\) is P — J2f=i ^^(QT)^ where A > 0, and , i = 1, . . . ,K, 
are the corresponding optimal covariance matrices for a fixed A in (27). 

Proof: The proof is provided in Appendix D. ■ 

Lemma 2 indicates that the value of A should increase, if ^^j^ ^^{QT) > ^^'^ "^^^^ versa. 
We are now ready to present our algorithm for solving Problem 4. 
Decoupled Iterative Power Allocation (DIPA) Algorithm : 

1) Initialize Amin and A^ax; 

2) repeat 

a) A = (Amin + Aniax)/2 

b) repeat, initialize Q^(0), ■ • • , Q^(0), n = 1 

for i = 1, ■ ■ ■ , K 



end for 

n = n + 1, 



QT{n-l)+tV^^lL 



c) until Q™ for k = 1, ■ ■ ■ , K converge, i.e., || VqLl|P < e for a small preset e. 

d) if Eti tr(Qr) > then A^in = A, elseif E^=i HQT) < then A^ax = A; 
3) until I Amin - AmaxI < e, 

where e > is a constant. The following proposition shows the convergence property of the 
DIPA algorithm. 

Proposition 4: The DIPA algorithm converges to an optimal set of the MAC transmit signal 
covariance matrices. 

Proof: The DIPA algorithm consists of the inner and outer loops. The inner loop is to 
compute Q™ for z = 1, ■ ■ ■ , A'. In each iterative step of the inner loop, we update by fixing 
other Q™ with j ^ i, and compute the corresponding gradient. The inner loop uses the gradient 
ascent algorithm, which converges to the optimal value due to its nondecreasing property and 
the convexity of the objective function. The outer loop is to compute the optimal Lagrangian 
multiplier A in (28). Due to the convexity of the dual objective function [23], there is a unique 
A achieving the optimal solution in (28). Hence, we can use an efficient one dimensional line 
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bisection search ( [19], [18]). ■ 
Remark 3: In the previous work on the sum rate maximization [19] [17] [18], the covariance 
matrix of each user is the same as the single user water-filling covariance matrix in a point-to- 
point link with multiuser interference being treated as noise [24] . However, for the weighted sum 
rate maximization problem, the optimal solution does not possess a water-filling structure. Thus, 
our DIPA algorithm does not obey the water-filling principle. In Section VI, Example 1 compares 
the water-filling algorithm with the DIPA algorithm. Notably, the formulation of Problem 4 is 
similar to the weighted sum rate problem for the dual MIMO MAC in [9]. The algorithm proposed 
therein to handle the dual MIMO MAC problem is based on gradient projection method [21]. 
The difference between our DIPA algorithm and the algorithm in [9] is just like the difference 
between the algorithms in [17] and [18]. 

The DIPA algorithm is an efficient algorithm to obtain the optimal transmit covariance matrix 
of the dual MIMO MAC (Problem 4). Moreover, the optimal solution to Problem 3 can be 
obtained via the MAC-to-BC covariance matrix mapping algorithm presented in the next section. 

IV. MAC-TO-BC Covariance Matrix Mapping 

A covariance matrix mapping algorithm was developed in [12]. However, this algorithm works 
for the sum rate maximization problem under a single sum power constraint, and is not applicable 
to a weighted sum rate problem under multiple constraints. In the following, we develop a 
covariance matrix mapping algorithm, which computes the BC covariance matrices via the 
dual MAC covariance matrices such that two channels yield a same weighted sum rate. 

In the MIMO-MAC, according to (11), the transmit beamforming vectors Vij can be obtained 
by the eigenvalue decomposition. The corresponding receive beamforming vector at the BS, Uij, 
is obtained by using the minimum mean square error (MMSE) algorithm: 

i-l Nr j-1 

= X] ^k,iHlvk,ivl iHk + q^H\vi^iv\^Hi + H\vi^j, (31) 

k=l 1=1 1=1 

where a is a normalized factor such that = 1. Throughout the proof of Proposition 3, we 

can see that when the same optimal solutions are achieved the primal BC and the dual MAC share 
the same beamforming vectors Uij and Vij. Hence, the transmit beamforming vectors of the BC 
are just the receive beamforming vectors of the dual MAC, and the receive beamforming vectors 
of the BC are the transmit beamforming vectors of the dual MAC. Thus, to obtain the transmit 
signal covariance matrix of SUj for the BC, we only need to compute the power allocated to each 
data stream. Due to Corollary 1, the dual MAC and the BC can achieve the same SINR region. 
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i.e., SINR^^ = SINRj" . Thus, for the BC, the power allocated to the beamforming direction Ui^j 
can be obtained by 

SINK- (eL+1 E^iPkA^Mv^.f' + E£+iP.d<.^^^..P + ^0 

For the BC, the encoding order is the reverse of the decoding order of the MAC. Thus, pK,Nr 
is computed first, pK,Nr-i is computed second, and so on, in the decreasing order of the data 
stream index and the user index. 

After computing the power for all the beamforming vectors, we obtain the signal covari- 
ance matrix from the BS to SUj, Q\ = ^^=iPi,jUiju\ ■. The aforedescribed process can be 
summarized as the following algorithm. 
MAC-to-BC Covariance Matrix Mapping Algorithm: 

1) Compute and Vi^j through eigenvalue decomposition: Q™ = VjAj V| = Yl!j=i ^ij^ij^l/^ 

2) Use the MMSE algorithm to obtain the optimal receiver beamforming vector Uij and 

sink;;;.; 

3) Compute pij through (32) according to the duality between the BC and the MAC; 

4) Compute Q- = Ylf=iPi,j'f^i,M,j- 

It should be noted that even though an explicit algorithm is not given, the paper [20] has 
mentioned the idea behind the above algorithm. The MAC-to-BC covariance matrix mapping 
allows us to obtain the optimal BC covariance matrices for Problem 3 by solving Problem 4. 

V. A Complete Solution to the CR MIMO-BC Weighted Sum Rate Problem 

We are now ready to present a complete algorithm to solve Problem 2. The Lagrangian dual 
objective function of Problem 2 can be rewritten as follows 

K 

9{(lt,qu)= max y^Wir^, (33) 
{Q^}fli:QVo 

where the maximization is subject to the constraint qt ( J2f=i ^'^oQ^^o — Pt) + Qu ( J2f=i ^^(Qi ) — 
Pu) < 0. Problem 2 is equivalent to the following problem 

min g{qt, g„), subject to qt > and qu > 0. 

<lt,qu 

Applying the BC-MAC duality in Section II-B and the DIPA algorithm in Section III, g{qt, qu) 
can be obtained. The remaining task is to determine the optimal qt and g^. Since g{qt, qu) is not 
necessarily differentiable, we search the optimal qt and g„ through the subgradient algorithm; 
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that is, in each iterative step, we update the vector [qt, g„] according to the subgradient direction 
s = [si,S2] of giquQu)- 

Lemma 3: The subgradient of giquQu) is [Pt - Yu'Li^lQ^ho, Pu - Ylf=i^AQ\)\^ where 
Qt > 0,qu > 0, and , i = 1, . . . , K, are the corresponding optimal covariance matrices for 
the problem (33). 

Proof: The proof is given in Appendix E. ■ 
It has been shown in [25] that with a constant step size, the subgradient algorithm converges 
to a value that is within a small range of the optimal value, i.e., 

lim - q:\ < e, and, lim |gW - qll < e, (34) 

n— »oo n—>QO 

(n) (n) 

where q^ and g* denote the optimal values, and ql and qu denote the values of qt and g„ at 
the nth step of the subgradient algorithm, respectively. This implies that the subgradient method 
finds an e-suboptimal point within a finite number of steps. The number e is a decreasing function 
of the step size. Moreover, if the diminishing step size rule, e.g., the square summable but not 
summable step size, is applied, the algorithm is guaranteed to converge to the optimal value. 

We next describe the algorithm to solve Problem 2 as follows. 
Subgradient Iterative Power Allocation (SIPA) Algorithm : 

1) Initialization: qt^\ qu \ n = 1, 

2) repeat 

2a) Find the optimal solution of the dual MAC Problem 4 through the DIPA algorithm; 
2b) Find the solution of the BC problem (33) through the MAC-to-BC mapping algorithm; 
2c) Update g^"^ and gi"^ through a subgradient algorithm g^"^^^ = g^"^ f^lQ^^o — 

Pt),q^u^'^=qt^+tiEtMQ')-Pu), 
2d) n = n + l 

3) Stop when |gJ"^(Eii /^IQ'^o - Pt)\ < e and |gi"^(E£i tr(Q- ) - Pu)\ < e are satisfied 
simultaneously, 

where t denotes the step size of the subgradient algorithm. As a summary, the flow chart of 
the SIPA algorithm is depicted in Fig. 3. We shows that the SIPA algorithm converges to the 
optimal solution of Problem 1 in the following proposition. 

Proposition 5: The SIPA algorithm converges to the globally optimal solution of Problem 1 . 
Proof: The Lagrangian function of Problem 1 is given by 

K K K 

i=l i=l i=l 
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and the Lagrangian function of Problem 2 is given by 

K K K 

U{Q\,...,Q)„X,qt,qu)=J2 " ^ ( E ^oQ'ho - Pt) - qu{J2 ^<Q') " ) • (36) 



i=l i=l 1=1 



Let qt, qu. A, and Qj be the optimal values of Li{Q^, . . . , Q^, A, gj, g„), when the algorithm 
converges. We thus have 

9Li(Q5,...,Q^,A,gj,g„) 



= 0, 



kt(I]£i ^iQi^o - Pt)\ = 0, and |g«(E£itr(Qi) - Pu)\ = 0. This means that is a locally 
optimal solution. 

According to (35), if we select Ai = Xqt, A2 = Ag«, and = Q^, then Ai, A2, and satisfy 
the KKT conditions of Problem 1 and thus are the locally optimal variables. 

Suppose that there exists an optimal set of Ai, A2, and such that L{Q^, . . . , Qk, Ai, A2) > 
L{Q^, . . . , Qj^, Ai, A2). Clearly, this optimal set of Ai, A2, and satisfy the KKT conditions 
of Problem 1 . In the sequel, we will derive a contradiction. 

First, we can write 

L(Qi, ■ ■ ■ , Qi^, Ai, A2) > ■ ■ ■ , Q^, Ai, A2). (37) 

Suppose that (37) does not hold, i.e., L{Qi, ■ ■ ■ , Q^, Ai, A2) < L{Qi, ■ ■ ■ , Qj^, Ai, A2). Then, 
according to the BC-MAC duality in Section II-B, an objective value of (7) which is larger than 
L{Qi, ■ ■ ■ , Qk^ Ai, A2), can be found for the fixed qt and g„. However, from Proposition 4, the 
DIPA algorithm converges the optimal solution. It is a contradiction. 

Secondly, according to the KKT conditions of Problem 1, we have 

K 

Xi{J2hiQ\-Pt)=0, (38) 

i=l 
K 

A2(5^tr(Q.)-Pu) =0. (39) 

i=l 

We thus can write: 

L(Qi, ■■■,Qk, Ai, A2) > ■ ■ ■ , Q^, Xu k). (40) 

Combining (40) and (37), we have 

L(Qi, ■ ■ ■ , Qi^, Ai, A2) > L{Q,, ■ ■ ■ , Qk, Ai, A2). (41) 
This contradicts with our previous assumption. ■ 
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Remark 4: The algorithm can be extended to the multiple PU case in the following manner. 
Assume that there are N PUs. Problem 2 becomes 

K 

min max y'wir^, 

N K K ^ ' 

subject to ^i^i ( H K,Q>o, - Pt,j) +qu{Y, • ) " Pu) < 0, 

j=l 1=1 i=l 

where qtj is the auxiliary variable for the jth PU, hoj is the channel response from the BS to 
the jth PU, and Ptj is the interference threshold of the jth PU. The role of auxiliary variables 
qtj is similar to that of qt in the single PU case. It is thus straightforward to modify the SIPA 
algorithm to solve the problem for the multiple PU case. Moreover, it should be noted that the 
multiple interference constraints of the problem (42) can be transformed to the per-antenna power 
constraints [20] by setting ho,j, j = 1, ■ ■ ■ , Nt, to be the jth column of the identity matrix. Not 
limited by the sum rate maximization problem with interference power constraints, the method 
proposed in this paper can be easily applied to solve the transmitter optimization problem (e.g. 
beamforming optimization) for MIMO BC with multiple arbitrary linear power constraints. 

VL Simulation Results 

In this section, we provide the simulation results to show the effectiveness of the proposed 
algorithm. In the simulations, for simplicity, we assume that the BS is at the same distance, li, to 
all SUs, and the same distance, /2"', to PU„. In the single PU case, we will drop the superscript 
and simply use notation I2. Suppose that the same path loss model can be used to describe the 
transmissions from the BS to the SUs and to the PUs, and the pass loss exponent is 4. The 
elements of matrix H are assumed to be circularly symmetric complex Gaussian (CSCG) RVs 
with mean zero and variance 1, and ho can be modeled as ho = {h/hYcin, where a„ is a A^t x 1 
vector whose elements are CSCG RVs with mean zero and variance 1. The noise covariance 
matrix at the BS is assumed to be the identity matrix, and the sum power and interference power 
are defined in dB relative to the noise power, and Pt is chosen to be dB. For all cases, we 
choose li = I2, except for explicitly stated. 

Example 1: In Fig. 4, we examine the validity of the DIPA algorithm. In this example, we 
choose K = 1 (a single SU case), Nt = A, = 4, and Pu = 10 dB. It is well known that 
the optimal transmit signal covariance matrix can be obtained through the water-filling principle 
[24]. As can be observed from Fig. 4, in several iterations, the DIPA algorithm converges to the 
optimal solution obtained by using the water-filling principle. 
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Example 2: In Fig. 5, we show the convergence property of the DIPA algorithm. In this 
example, we choose K = 20 and P„ = 10 dB. It can be observed from this figure that the 
algorithm converges to the optimal solution within several iteration steps. 

Example 3: In Figs. 6 and 7, we consider a SU MIMO-BC network with K = 5, Nt = 5, 
Nr = 3, and P„ = 13 dB. In this example, the SUs with wi = 5 and Wi = 1, i = 2, . . . , K 
are assumed to share the same spectrum band with two PUs. Fig. 6 plots the weighted sum rate 
versus the number of iterations of the SIPA algorithm for step sizes t = 0.1 and t = 0.01. As 
can be seen from the figure, the step size affects the accuracy and convergence speed of the 
algorithm. Fig. 7 plots the sum power at the BS and the interference power at the PUs versus 
the number of iterations. It can be seen from the figure that the sum power and the interference 
power approach to = 13 dB and P^ = dB respectively when the SIPA algorithm converges. 
This implies that the sum power and interference constraints are satisfied with equalities when 
the SIPA algorithm converges. 

Example 4: Fig. 8 plots the achievable sum rates versus the sum power in the single PU case 
and the case with no PU. We choose K = 5, Nt = 5, and N^- = 3. As can be seen from Fig. 
8, in the low sum power regime, the achievable sum rate in the case with no PU is quite close 
to the one in the single PU case while in the high sum power regime, the achievable sum rate 
in the case with no PU is much higher than the one in the single PU case. This is because the 
additional constraint reduces the degrees of freedom of the system. 

Example 5: In this example, we consider the influence of the interference constraint on the 
achievable sum rate of the SUs. In this example, Nt = 5, K = 5, and Nr = 3. The sum power 
constraint for the BS is assumed to be 15 dB and 20 dB. Fig. 9 compares the sum rate achieved 
in a PU case with one achieved in the case with no PU as I2/I1 varies from 1 to 12. It can be 
observed from the figure that the achievable sum rate increases as the PU moves away from the 
BS, and the influence of the PU reduces to zero after the I2/I1 is larger than a certain threshold. 

VII. Conclusions 

In this paper, we developed a new BC-MAC duality result, which can be viewed as an extension 
of existing dual results developed under either a sum power constraint or per-antenna power 
constraints. Exploiting this duality result, we proposed an efficient algorithm to solve the CR 
MIMO-BC weighted sum rate maximization problem. We further showed that the proposed 
algorithm converges to the globally optimal solution. 
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Appendix 



A. Lemma 4 and its proof: The following lemma describes an important property that will be 
used in the proof of other lemmas. 

Lemma 4: For fixed qt and g„, the maximum weighted sum rate in (5) is achieved when the 
constraint (6) is satisfied with equality. 

Proof: We here adopt the DPC scheme, which is a capacity achieving strategy for the 
MIMO-BC [16]. Let the permutation vr represent the encoding order when the optimal solution 
is achieved. Assume that SU,r(i) is encoded first such that the signal of SU^(i) is noncausally 
known to the BS before the signals from the other SUs are encoded. Thus, in the DPC scheme 
the signal from SU7r(i) has no impact on the rates achieved by the other SUs. We prove this 
lemma by contradiction. 

Suppose that is the optimal signal covariance matrix of SU7r(i). Assume that the con- 

straint (6) is satisfied with a strict inequality when the optimal solution is achieved. Thus, we 
can always find an e > such that 

K K 



i=2 



i=2 



K 



K 



i=l 



i=l 



Moreover, the rate achieved by user 7r(l) in the MIMO-BC can be written as 



log 



Due to the positive semi-definiteness property of Q\, we have 

K 

log 



log 
log 



1=2 
K 



i=2 



K 



K 



/2 



i=2 



i=2 



April 22, 2008 



DRAFT 



21 



K 



= log 
= log 
>log 
= log 



1=2 
K 

i=2 
K 

i=2 
K 



+ log 
+ log 
+ log 



I + G\Ql^,^ + eI)G 



/ + A + eS 



/ + A 



(44) 



4 = 1 



(45) 



where G = + J2f=2-^'^wQn{i)Hi{i))~^^'^^ ^^'^ ^ ^^'^ ^ diagonal matrices. Eq. 

(44) is due to the fact that the optimal covariance matrix for a MIMO has the water-filling 
structure [19] [24], i.e., if we apply singular value decomposition to G, G = VSU, where V 
and U are unitary matrices, and S is a diagonal matrix, then the optimal Q^(^i) can be written 
as = U^RU, where i? is a diagonal matrix. Thus, we have A = SRS and S = 55. 

According to (45) and (43), + eJ is a better solution for the sum rate problem than 

Q^(^), which contradicts with the assumption. Therefore, the constraint must be satisfied with 
equality. ■ 

B. Proof of Proposition 1: The proof consists of two parts. In the first part, we show that either 
optimal solution is feasible for both problems. In the second part, we show that Problem 1 and 
Problem 2 have the same solution. 

The Lagrangian function of Problem 1 is 

K K K 

LiiQ\,---, Qk, a*, A„) = J2 - A, ( ^ hlQ'lK - Pt) - \u ( ^ tr(Q^) - P„) , (46) 



i=l 



i=l 



i=l 



where At and A„ are the Lagrangian multipliers. The optimal objective value is 

min max Li{Q\, ■ ■ ■ , Q)^, \t, \u). 



(47) 



Assume the optimal variables are Xt, A« and Q\, - ■ ■ , Q^, and the corresponding optimal value 
is C. 

The Lagrangian function of Problem 2 is: 

K K 

L2{Q\r--. Qk, qt, Qu, A) = w,r^ - A (g* ( ^IQ ■ -Pt)+qu{Yl tr(Q^) - P„) ) , (48) 

i i=l 1=1 
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where A is the Lagrangian muhiplier. The optimal objective value is 



min max LiiQ^, ■ ■ ■ , Qk, Qt, Qu, >^)- (49) 



Suppose that the optimal variables are qt, (ju, A, and , i = 1, . . . , i^, and the corresponding 
optimal objective value is C. We just need to prove C = C. 

We now present the first part of the proof. According to the KKT condition of Problem 2, we 
have 

dUKQl-.^.^Q^AtAu) ^tgj^^ _ p^) = 0, (50) 

^* 1=1 

?MMk_^k^.X(f:.r(Qb-/'.)=0. (51) 

Recall that the Lagrangian multiplier A is non-negative. Furthermore, if A = 0, we have qt{J2f=i 
Pt)+Qui^f=i tr(Oj) < from the KKT conditions. This contradicts with Lemma 4. Thus, 
we always have A > and can readily conclude that Yld=i ^\Q\^o = Pt and Yl!i=i ^AQ\) = Pu 
are satisfied simultaneously. The optimal solution of Problem 2 is also a feasible solution of 
Problem 1 . On the other hand, it is obvious that the feasible solution for Problem 1 is also the 
feasible solution for Problem 2. 

We next prove the second part by using contradiction. Let us first suppose C > C . For (48), 
if we select Q\ = q\ fox i = 1, . . . , i^, A = 1, qt = \t and g„ = A„, then L2 = C > C. li 
contradicts to the fact that C is the optimal objective value for (49). 

We now assume C <C. Recall that A 7^ 0, for (48). If we select Q\ = q\ fox i = I, . . . , K , 
Xt = Xqt and A„ = Ag„, then Li = C > C, which contradicts with the fact that C is the optimal 
objective value for (47). 

Therefore, the optimal solutions for Problem 2 and Problem 1 are the same. ■ 

C. Proof of Lemma 1: According to previous discussions, the signal from each SU is divided 
into several data streams. We now show that the optimal encoding order of these data streams 
are arbitrary. It is well known that the optimal objective value of the MAC equally weighted sum 
rate problem can be achieved by adopting any ordering [19] [17] [18]; that is, when all the users 
have the same weights, the optimal solution of the weighted sum rate maximization problem is 
independent of the decoding order. Analogously, the data streams within a SU share the same 
weight. Thus, an arbitrary encoding order of those data streams within a SU can achieve the 
optimal solution. ■ 
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D. Proof of Lemma 2: Let s be the sub-gradient of g{\). For a given A > 0, the subgradient s 
of giX) satisfies g{\) > g{\) + s(A — A), where A is any feasible value. Let Q™, i = 1, . . . , J^T, 
be the optimal covariance matrices in (27) for A = A, and « = 1, . . . , -R', be the optimal 
covariance matrices in (27) for A = A. We express g{\) as 

K 

g{\) = max (/(Q-, ■ ■ ■ , Q?^) - A(5^tr(Q-) - P) 

= /(Qr,---,Qx)-A(X:tr(Or)-^ 

j=i 

>/(Qr,---,Q^)-A(x:tr(Qr)-p) 

i=l 

= /(or, ■ ■ ■ , - A ( E tr(Qr)-i^)+ A ( x: tr(or) -^) - a ( f: ^or) - 

i=l i=l i=l 

K 

= g{\)+{P-Y.tr{Q:))C\-\), 

i=l 

where s := P — J2f=i ^^(QT) the subgradient of g{X). This concludes the proof. ■ 

E. Proof of Lemma 3: The subgradient s of g{qt, qu) satisfies g{qt, qu) > g{qt, Qu) + (fe, Qu] - 
[qt, qy\) ■ , where [g^, g„] is any feasible vector. Let z = 1, . . . , /T, be the optimal matrices 
of the problem (33) for qt = qt and = g^, and let i = 1, . . . , be the optimal matrices 
of the problem (33) for qt = qt and g„ = g„. We express g{qt., qu) as 

M 

9{Qt, Qu) = max ^ ^ u/^r,^ (52) 
Ql'"' 'Qk i=i 

M K _ K _ 



i=l i=l i=l 

M K K 

> ^.^t - A (g* ( 5^ hlQX - Pt) +qu{Y. • ) - Pu) ) (54) 

i=l i=l i=l 

M K K 

y^^r'^ - A (qt ( ^''oQX - Pt) +qu{Y '<Q') " ) 

i=l i=l i=l 

K K K K 

+ Afe [Y hlQX-Pt)+qu iY ^<Q')-Pu)) - Afe iY hlQX-Pt)+<lu iY tr(Q • )-^«)) 



1=1 1=1 1=1 i=l 
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K K 



--9{qu Qu) + [Y1 '^IQ'^^'^o - Pt) {Xqt - Mt) + ( tr(Q') - P„) (Ag„ - Ag„) 

1=1 i=l 
K 

--9{Qu Qu) + (^Y hlQ^ho - Pt) {Xqt - Xqt + Mt - Mt) 

i=l 

K 

+ ( 5Z ^^(^i ) ~ (^^» ~ -^^^ + ~ '^^") 

K K 

-9{qt, Qu) + ( ^^oQ\ho - Pt) {Xqt - Xqt) + ( ^ ^oQX - Pt) {Xqt - Xq, 

i=l i=l 

K 

+ ( Y - P^ (Xqu - Xqu) + ( Y - - ^^-) 



i=l 



=9{qt. qu) + ( 5^ hiQ% - P^) {Xqt - Xqt) + ( ^ tr(Q') - P„) (Ag„ - Ag„) (55) 

i=l i=l 

=9{(it, qu) + X{[qt, qu] - [qt, qu]) ■ s^, 

where s := [Pt — Yl,f=i^\Q\^o-, Pu — Ylf=i^AQl)]- Eq- (53) is due to the fact that the dual 
objective function of the problem (33), and r°. A, and Q\ are the optimal variables for the fixed 
qt and g„. The inequality (54) is because , z = 1, . . . , are the optimal signal covariance 
matrices for the fixed qt and g„. The equality (55) is due to Lemma 4. Thus, s is the subgradient 
offi'(gt,g«). ■ 
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Fig. 1. The system model for MIMO-BC based cognitive radio networks. Tliere are K SUs. The BS of the SUs has Nt 
transmit antennas, and each SU is equipped with Nr receive antennas. 
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Fig. 2. The system models for Problem 3 and Problem 4, where qt and q„ are constant, and Ro = hohl. 
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Fig. 3. The flow chart for the SIPA algorithm, where and Q" denote the transmit signal covariance matrices of SUi 

for the BC and MAC at the nth step, respectively. 
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Fig. 4. Comparison of the optimal achievable rates obtained by the DIPA and the water-filling algorithm in a MIMO channel 
(Nt = Nr = 4, K ^ 1 and P„=10 dB). 
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Iteration step 

Fig. 5. Convergence behavior of the DIPA algorithm (K = 20 and P„ — 10 dB). 
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Fig. 6. Convergence behavior of the SIPA algorithm {Nt = 5, K = 5, Nr — 3, wi = 5, and Wi — 1, for i ^ 1). 
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Fig. 7. The convergence behavior of the sum power at the BS and the interference at the PU for the SIPA algorithm (Nt = 5, 

K = 5, Nr = 3, wi — 5, and 1^ = 1 with i ^ 1). 
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Sum Power (dB) 

Fig. 8. Achievable sum rates versus sum power in the single PU case and the case with no PU (Nt = 5, K — 5, Nr = 3). 




Fig. 9. Achievable sum rates versus the ratio of h/h using the SIPA algorithm (Nt = 5, K = 5, Nr = 3). 
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